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Abstract 

Both ecological systems and convective fluid systems are examples 
of open systems which operate far-from-equilibrium. This article demon- 
strates that there is a correspondence between a resource-consumer chemo- 
stat ecosystem and the Rayleigh-Benard (RB) convective fluid system. 
The Lorenz dynamics of the RB system can be translated into an ecosys- 
tem dynamics. Not only is there a correspondence between the dynami- 
cal equations, also the physical interpretations show interesting analogies. 
By using this fluid-ecosystem analogy, we are able to derive the correct 
value of the size of convection rolls by competitive fltness arguments bor- 
rowed from ecology. We finally conjecture that the Lorenz dynamics can 
be extended to describe more complex convection patterns that resemble 
ecological predation. 

Keywords: Rayleigh-Benard system; Lorenz model; resource-consumer chemo- 
stat; ecosystem metabolism; thermodynamics 

1 Introduction 

In a Rayleigh-Benard experiment, a horizontal viscous fluid layer is heated from 
below. When the temperature difference between upper and lower sides is small, 
heat transfer solely occurs through thermal conduction. Yet once beyond a 
critical temperature difference, a regular pattern of convection cells or rolls 
emerges (Benard, 1901). This sudden shift from conduction to convection is 
referred to as the Rayleigh-Benard (RB) instability, and is often quoted as an 
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archetypal example of self-organization in non-equilibrium systems (Nicolis and 
Prigogine, 1989; Prigogine, 1967). 

Intuitively, it makes sense to try to apply the concept of self-organization in 
non-equilibrium systems to ecologicaly systems, as there are similarities between 
ecological and physical systems. Like the Rayleigh-Benard set-up, ecological 
systems are open systems that receive a throughput of energy and/or mass 
via coupling to two environments (Morowitz, 1968; Schrodinger, 1944). These 
two environments are typically large reservoirs and they drive the system from 
equilibrium. Consider the example of a laboratory chemostat ecosystem (Smith 
and Waltman, 1995). This is a prime example of a chemotrophic ecosystem 
whereby a resource of energetic high quality chemical substrate is pumped from 
a reservoir into the system. In the ecosystem this resource is degraded into low 
quality waste products which are emitted to the waste reservoir. When there 
is low feeding of resource, no biota can survive, and the resource is degraded 
by abiotic processes only. But when the feeding is above a critical threshold, 
biota can survive by consuming the resource. There is a sudden shift from a 
lifeless to a living statfH- In other words, the energetic quality difference between 
incoming and outgoing chemical substrates is exploited by various abiotic and 
biotic processes. The latter biotic processes contain the biomass synthesis and 
turnover of consumer micro-organisms feeding on the resource. 

So it is tempting to look for a deeper connection. Can one compare bi- 
ological processing with convection? Both mechanisms involve self-organizing 
structures, biological cells or convection cells, that can only survive after a crit- 
ical threshold. Both energetic pathways, biotic resource conversion and thermal 
convection, degrade energy from high quality to low quality form. And these en- 
ergetic pathways are additional to the abiotic conversion or thermal conduction 
processes of the background. 

Here, our ambition is to examine the link between ecological processes and 
convective fluid motions in a quantitative way. The first part of this article con- 
tains a highly intriguing result: The mathematical expressions of the resource- 
consumer chemostat ecosystem dynamics are exactly the same as the dynamics 
that describes the basics of the Rayleigh-Benard system. Furthermore, not only 
are the mathematical equations identical, also the physical/ecological interpre- 
tations give appealing results. Particularly, by looking at the energetic pathways 
of the ecosystem, the ecological quantities can be mapped to the quantities used 
in the fluid system and vice versa. 

The second part tries to extend the correspondence between fluid convection 
and ecosystem functioning to include new processes. We will study two exten- 
sions. First, one can look at ecological competition and translate the notion 
of competitive fitness to the fluid system. The convection cells are in 'Dar- 
winian competition' with each other and the fittest ones will survive. One can 
generalize the Lorenz model to include this fluid competition. As the size of a 
convection cell will depend on the fitness measure, we will demonstrate that the 
mathematical identity of the ecological and the fluid dynamics predicts the ex- 
perimentally correct size of the cells at the onset of convection. Second, one can 
look at ecological predation. Translating this notion to the fluid system leads to 
a new conjecture to extend the Lorenz model in order to describe more complex 

^Strictly speaking, it is rather a distinction between guaranteed extinction and survival. 
We study what will happen with an organism which is released in the ecosystem. Our results 
should not be interpreted as the solution for the origin of life. 
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convection patterns. These new patterns only appear when the system is driven 
beyond a second critical value for the temperature gradient. The 'predatory 
behavior' in the convective fluid system leads us to a conjecture which we will 
not prove, but will be successfully tested by looking at the energy dissipation. 

2 The Rayleigh-Benard convection system 

Let us start by deriving the dynamics that describes the Rayleigh-Benard (RB) 
convective fluid system, named after Benard (190f) and Rayleigh (1916) who 
were the first to study this system experimentally and theoretically. A full 
mathematical treatment of thermal convection requires the combined solution 
of the heat transport, Navier-Stokes and incompressibility equations, resulting 
in a set of five coupled non-linear partial differential equations (Chandrasekhar, 
1961; Rayleigh, 1916). 

Rather than solving this full set, we employ the approximation adopted by 
Lorenz (1963), which became famous as it gave an impulse to the development 
of chaos theory. The model describes the lowest modes of an expansion of 
the temperature and velocity fields for a RB system with free-free boundary 
conditions (see e.g. Getling 1998). In appendix iBl the derivation of the Lorenz 
system is given in a way that will suit our further discussion. A non-linear set 
of three ordinary differential equations is obtained, with three variables (see 
fig. [T]): X measures the rotational rate of the rolls and represents the maximal 
velocity at the bottom of the rolls. Y and Z are temperature deviations, where 
the linear profile of the conduction state is taken as a reference. 




X 



Figure 1: The profiles and the variables X, Y and Z. The temperature (T) 
and horizontal velocity (vx) profiles at three vertical sections (dashed lines) are 
shown. These vertical sections are parallel with the axes of the convection rolls, 
where the fluid is moving up, moving horizontal or moving down. The thin 
linear profiles correspond with the conduction state, the thick profiles with the 
convection state. As indicated, Y and Z are temperature deviations and X is 
the velocity at the bottom of a roll. 

With these three variables, the XYZ Lorenz system is rich enough to describe 
the Rayleigh-Benard instability, the sudden shift from conduction to convection. 
But there is an even simpler model, the XZ system with only two variables, that 
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is rich enough as well. It is this XZ model that allows us to make the corre- 
spondence. Roughly speaking, we will perform an averaging over the horizontal 
directions, such that only the average vertical profile remains. As Y is the tem- 
perature deviation in horizontal direction, it is this variable that will disappear 
after the averaging. Specifically, this is done by making the pseudo steady state 
assumption dY/dt — for the variable Y. The latter becomes a constant and 
the dynamics turns into: 

dX f a^h'^ga [a^ + 1)71"^ v\ 2a^hqa 



dt V(a^ + l)^'r^X ) (a2 + l)27rx 

dZ_ ^ a^Tr(3 ^ Q^^r^ ^n^X /„x 

dt 2(a2 + l)hx (a2 + l)h^x ' 

with h the height of the fluid layer, a a geometric factor such that h/a is the 
width of the straight convection rolls, a the thermal expansion coefficient, g the 
gravitational acceleration, x the heat conduction coefficient, v the kinematic 
viscosity and 

^=-^^ (3) 

the temperature gradient. This important quantity is the thermodynamic gra- 
dient that drives the system out of equilibrium. Th is the high temperature of 
the heat reservoir below the fluid layer and is the low temperature of the 
heat reservoir above the layer. 

For further reference, we will also need measures for the temperatures at the 
middle and the lower halve of the fluid layer. Define 



^ (4) 



T^ + T? 
"2 

as the horizontally average temperature at height h/2, and 



T„ . a±Z|L^, (5, 

Fig. [2] shows the interpretation of Th as a temperature measure for a linearized 
temperature profile in the lower half of the fluid layer. (Due to symmetry in 
the approximation leading to the Lorenz system, we will not have to include the 
upper half of the fluid layer.) 



3 The resource-consumer ecosystem 

Next, we discuss the ecosystem model, which is in essence a simple chemotrophic 
resource-consumer food web model, one of the mainstay models of ecology (e.g. 
Yodzis and Innes, 1992). Consumer organisms are feeding on some food re- 
source (R), which is partly converted to consumer biomass (C) and partly to 
waste product (W). For reference, one can think of a chemostat set-up where 
a chemical reactor tank contains a monoculture of micro-organisms that are 
feeding on a chemical substrate like methane or glucose, while respiring C02- 

Figure [3] gives a schematic overview of the ecosystem coupled with the two 
environments (denoted with the superscripts 0). There are two environmental 
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Figure 2: The definition of Th- The lower half of the fluid layer is shown, with 
the vertical (horizontally averaged) temperature profiles in the conduction state 
(thin line), the convection state (thicker line) and the 'linearized' convection 
profile (thickest line). The variables Z and Th are at height h/A. 




Figure 3: The ecosystem fiow chart. The resource-consumer-waste ecosystem 
coupled with the two environments. The different fiuxes are discussed in ap- 
pendix]^ The color denotes the 'energetic quality' of the substances, from high 
(red) to low (fight yellow). 



compartments, the resource at constant concentration and the waste at 
constant concentration Cy^) and three ecosystem compartments, with variable 
concentrations Cr, Cc and Cw for the resource, the consumer biomass and the 
waste respectively. In appendix |^ the complete dynamics of the resource-con- 
sumer-waste (RCW) ecosystem is given, explaining the fiuxes F between the 
compartments. 

However, as we will see, the correspondence only works in a limiting case, 
whereby roughly speaking we will average over the waste concentrations of the 
system and the environment. Specifically, this can be done by taking a very 
small relaxation time for the exchange of the waste between the ecosystem 
and the reservoir. This means that by studying the ecosystem at longer time 
scales than this relaxation time, the dynamics for W is forced to be in a pseudo 
steady state condition. Hence, W is no longer a variable and we end up with the 
resource-consumer (RC) model, with two dynamical equations for two variables: 

aR{Cl - Cr) - {kac + gcRCc){CR - C/i^e,), (6) 
qCRgCR{CR - C^/Keq)Cc - dcCc, (7) 



dt 
d 

di 



c 
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with, an the resource exchange rate parameter, kac the abiotic conversion (from 
R to W) rate parameter, gcR the consumer growth rate parameter, qcR the yield 
factor for the consumer growth, dc the consumer decay (biomass turnover) 
rate parameter, and Kgq the equihbrium constant for the chemical reaction 
(oxidation) from R to W which always slowly proceeds at the background. 

This is the well-known chemostat dynamics (Smith and Waltman, 1995), 
which is extended in two ways: First, abiotic conversion is included in terms of 
chemical oxidation with parameters kac and C^/Keq. Second, instead of the 
classical dependence of the growth on the resource Cr, the growth is now made 
dependent on Cr — C^r/Keq. This is done for thermodynamic consistency: at 
chemical equilibrium, biomass synthesis should also cease. 
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Table 1: Two corresponding models with analogous mechanisms for gradient 
degradation 

Let us summarize. Table [T] shows the observation that there are two sys- 
tems with analogous mechanisms for the degradation of a gradient, i.e. the 
transformation of high quality energy to low quality energy. The unstructured 
processes are the ground level mechanisms: abiotic conversion from resource to 
waste or thermal conduction from high temperature to low temperature. But 
above a certain critical threshold, a self-organization mechanism adds second 
level processes: biotic conversion or thermal convection. 

To study these systems, we introduced two models, each with three variables: 
the XYZ Lorenz model (one velocity X and two temperatures Y and Z) and 
the RCW ecosystem model (one biotic consumer C and two abiotic molecules R 
and W). These systems have different behavior, as the XYZ model has chaotic 
solutions whereas they are absent in the RCW model. However, there is a 
hidden correspondence which we will clarify in the next section. We have to 
make a pseudo steady state condition (an averaging) of the 'abiotic' variables Y 
and W, leading to the XZ model and the RC model P?]) . These models 
have only two variables and hence they are the most simple models to study 
a non-trivial behavior, the transition from an 'abiotic' to a 'biotic' state. The 
XZ system does not have chaotic solutions anymore, so it is possible that it 
is mathematically equivalent with the RC model. The trick is to rewrite the 
variables and the parameters to demonstrate this equivalence. To give a first 
hint, the basic observation is that the variables should be related as 



(9) 



Note that the quantities on the right hand side have dimensions of temperature. 
In the next section, we will also relate the parameters and discuss the physical 
interpretations of this correspondence 



4 The correspondence 

So it is time to write the dictionnary of the correspondence. Our final result is 
shown in table [5] at the end of this article. In order to reach our goal, we need to 
be able to consistently translate quantities from one system to the other. The 
redefinitions explained below enable us to write the simplified Lorenz dynamics 
as the ecosystem dynamics. 

First we will state the relation between the basic quantities, the concentra- 
tions and the temperatures, which is simply: 

c° ^ n, (10) 

Cb. ^ Th, (11) 
^ - K. (12) 



Keg 



These were derived by using ^ and the interpretation of Th ©. It explains 
why we can roughly interpret the resource as the heat energy. 
The consumer concentration is given by ([9]). As X is a velocity measure, is 
a measure for the kinetic energy of the convection rolls. This kinetic energy is 
consuming the heat energy resource. 

The yield and consumer growth parameters are written as 

[a^ + l)x 

As the gravitational field is causing the buoyancy force, this explains why g 
appears in qcR- Furthermore, these parameters depend on geometric factors, 
especially ft./a, the width of a convection roll. The importance of this depen- 
dence will be shown later. 

The abiotic exchange and abiotic conversion parameters are 

2^^'^ ^^^^ 

As these parameters are conduction coefficients, it is logical that they depend 
on the heat conduction coefficient x- In order that the analogy works, our 
ecosystems should have equal exchange and abiotic conversion parameter^ 



The reason is that we identified lllll l. and the distance between the lower side and height 
h/4 equals the distance from this height to the middle of the fluid layer. There is a possibility 
to have a more general correspondence, with an kaCj but then we will loose the relation 
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The final parameter is the biomass decay rate 



2{a' + ly^v 

dc ^ -f^ . (16) 

This explains why this decay is a kind of friction term. As mortality and viscous 
friction destroy the biological or convcctivc cells, a continuous feeding on the 
resource is required in order that these structures can survive. 

Having discussed the relations between variables and parameters of both 
systems, one can take a look at other ideas and concepts of one system and 
translate it to the other. A quantity that will become useful later is the ther- 
modynamic gradient that measures how far the system is out of equilibrium. It 
is given by the difference in energetic 'quality' of the two reservoirs. 

A° = C°-^-M = (tO_TO,). (17) 

The latter relation can be turned into a dimensionless measure, which is the well 
known Rayleigh number Ra in fluid systems. Another important dimensionless 
fluid quantity is the Prandtl number Pr. We can now see that they can be 
casted into their ecological analogs: 

Pr ^ "L^H . (19) 



X aR + KAC 



The geometric factors 



n = (21) 

will become important later on. 

This is the first part of our dictionarry. In the next section we will delve 

deeper into the physical analogies between both systems. In particular, we will 
look at the energy dissipation along the different energetic pathways. 



5 Energy flows along energetic pathways 

Our next challenge is to see whether the correspondence also works for the 
energy flows along the different pathways. Are the heat transport and the 
ecosystem metabolism connected? This question is not trivial, because even 
though the dynamical equations look the same, a priori it is not obvious that 
the thermodynamical expression for the heat transport is exactly the same term 
in the dynamical equations which corresponds with the ecosystem metabolism 
rate. Schneider and Kay (1994, fig 2a) used experimental data sets for the RB 
system to plot the total steady state vertical heat transport per unit horizontal 
area W* . (Steady states are denoted with a superscript *.) Our approach now 
allows us to write down a simple analytical expression for the heat transport in 
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the steady state, because it will be shown to be related with the total ecosystem 
metabolism (the total rate of waste production, see (|A-5[) ) 

Fem = {nAc + (1 - qcR)gcRCc){CR - Clr/K,,) + dcCc- (22) 

Our result will suit well the behavior as seen in the plot derived by Schneider 
and KajH. 

To calculate W* , observe that there is no energy accumulation in the fluid, 
and hence this heat transport is the same at every height. Therefore it equals 
the transport at height z — Q. At the bottom layer, the vertical fluid motion is 
zero, as is seen in the chosen boundary condition (|B-6p . Hence, at the bottom 
layer there is no vertical heat transport by fluid motion. The heat transport is 
given by the temperature gradient only, as for the conduction state. Taking a 
horizontal average, the Y-ieixn in the expansion (|B-lip drops out, leaving only 
the Z-term. This gives: 

rIT* 

{pocv)-^W* = -X^U=o (23) 
oz 

= x(/3+^) (24) 
= '^an{Cl-C*j,). (25) 

{po is the reference density and Cp is the heat capacity.) The latter expression 
gives the steady state resource exchange Fj^, which equals F^j^j (this is easily 
seen because there is no accumulation of ecosystem resource or biomass, and 
hence the net resourche exchange should equal the total conversion from resource 
to waste). 

In order to study the behavior of the ecosystem metabolism F^j^^ under 
different gradients, we need to solve the dynamics for the steady states. Scanning 
A*^ from zero to infinity, there is a critical value given by the bifurcation point 

^0 ^ (gfl + nAc)dc ^26) 
anqcRgcB^ 

For a value of A" < A°, we have only one stable steady state that is physically 
realistic (no negative concentrations) 

aR + KAC Keq 

Cc = 0. (28) 

Within this region, a stable population of consumers cannot be formed, and 
hence, only abiotic degradation takes place. However, if the resource input 
increases so that A" > Aj!, there is the possibility for the consumers to survive 
at a non-zero concentration. The above state becomes unstable, and the new 
stable solution becomes 

Q = -^ + ^, (29) 

qCRQCR J^eq 

^, _ qcR9CRaR^° - {aR + KAc)dc , . 

- d^R • ^^^^ 



^Schneider and Kay (1994) described a fluid layer with rigid-rigid boundary conditions. 
Therefore, our results can only be compared qualitatively, as our XZ model only works for 
systems with free-free boundary conditions. 
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Figure 4: The total steady state ecosystem metabolism F^j^j in function of 
the driving force A°, for specific parameter values. Red color denotes abiotic 
conversion only, the green line corresponds with biotic consumption. The dashed 
line corresponds with unstable states. 

Using these solutions, we can plot the steady state ecosystem metabolism 
^EM ^ function of the thermodynamic gradient A°, Figured The (qualita- 
tive) similarity with the Figure 2a in Schneider and Kay (1994) is obvious. The 
steady states which have only abiotic conversion are located at the so called 
thermodynamic branch, because this branch contains thermodynamic equilib- 
rium at zero gradient {F'^j^j = at A° = 0). In the RB system, these states 
correspond with thermal conduction. But above the bifurcation point, there is 
an exchange of stability: the thermodynamic branch states become unstable and 
new stable states arise. These are located at the so called dissipative branch, 
and they contain both abiotic and biotic degradation of resource. Once beyond 
the bifurcation, a viable consumer population can be established. Translated to 
the RB system, both conductive and convective heat transport processes appear 
and a viable 'kinetic energy population' is established. 

The above discussion shows the exact correspondence between two terms 
in the dynamics that describe the energy dissipation: the heat transport and 
the ecosystem metabolism. However, the argument was restricted to the steady 
state behavior. We will now give some other arguments to demonstrate that 
there is not only a formal mathematical equivalence of the RC and the XZ 
models, but that the terms in the dynamical equations correspond also physically 
with the different energetic pathways, Fig. [5] (compare with Fig. [3]). This 
correspondence of the energetic pathways of both systems is also valid in the 
transient states. 

First let us look at the exchange with the external reservoir (E). The fluid 
has a heat exchange with the heat reservoir at constant temperature T^. This 
exchange is due to heat conduction with coefficient x- The ecosystem has the 
same functioning: The variable Cn is in contact with the constant with 
exchange rate an, explaining the relation (fT5|) and ([25|) . 

Next, let us focus on the energetic pathways within the system. In our 
resource-consumer ecosystem, we have seen that there are basically three meta- 
bolic pathways for the consumption of the resource Cr (see Table [3] in appendix 
El). These are the three arrows arriving at the waste compartment C^/ in the 
figure. 

Also our fluid system has three equivalent heat transport and energy trans- 
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Figure 5: The correspondence between temperatures and concentrations, as well 
as the energetic pathways. The figure shows the lower half of the convective 
fluid layer, as in Fig. [21 with vertical energy transport. The boxes represent 
the heat energy compartments at three different heights, plus the kinetic energy 
compartment. In this way, the Rayleigh-Benard system is presented in terms of 
a simple resource-consumer food web. See text for more explanations. 

formation pathways (see e.g. (HH)) : 

• Ml: There is heat transport by conduction which is qualitatively given 
by Wcond oc x{Th — Tli)/{h/'i), and this is indeed proportional with the 
abiotic conversion 



• M2: There is direct heat transport by convection, i.e. heat energy from 
the lower reservoir is actively transported to the upper reservoir, without 
being turned in kinetic energy. This is a loss term for the transforma- 
tion of heat energy to kinetic energy. It is easily seen that this term is 
proportional with the consumer consumption (which is coupled with the 
consumer growth) 



because Cc is the kinetic energy. Using the dictionary, one can translate 
this expression into an analytical expression for the direct heat transport 
by convection Wdtr^conv 

• M3: There is heat production due to viscous dissipation of kinetic energy. 
This extra heat produced is also finally released in the cold temperature 
reservoir. It is the indirect heat transport by convection, as the heat energy 
is first turned into kinetic energy, and eventually released again as heat 
energy. In Kreuzer (1981), a derivation is given for this transformation 
rate of kinetic energy into heat energy: 



FAc^^iAciCR-C^/K,q). 



(31) 



(1 " qcR)FBc = (1 - qcR)gcR{CR - C\ 



(32) 
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which is indeed proportional with vX^ and hence with biomass decay 
Fed = dcCc (which in the steady state equals the consumer growth 
(IcrFbc)- 

To summarize, we have demonstrated a unique example of a correspon- 
dence between a biological and a physical system. The dynamical equations are 
equivalent and a dictionary was given between the different quantities. Also the 
physical interpretations (in terms of energetic pathways) of the different terms 
in the dynamical equations were proven to be analogous. This correspondence 
allowed us to calculate an analytical expression for the heat transport in the 
steady state of the RB system. As Fig. [5] shows, a simple resource-consumer 
food web arises in the fluid system. In the next two sections, we will take this 
analogy some steps further by expanding the fluid food web in two ways: First 
we will include competition at the first trophic level (the level of the consumers). 
Secondly, we will study longer food chains by including predation. In a sense, 
this approach allows us to use ecological concepts to extend the Lorenz dynam- 
ics in order to find new solutions (i.e. new convection patterns) for the fluid 
system. 



6 Competitive exclusion and fitness 

In ecology, there is the important idea that species can mutate and evolve, lead- 
ing to Darwinian competition between species. If we describe competition in our 
ecosystem by taking n different consumer species with growth rates gem, death 
rates dci and yields qcm, with i — 1, ...,7i, we can calculate the stable steady 
state and it appears that the species with the highest value of the competitive 
fitness 

r _ qcRigcRi 

da 

survives, the others go extinct. This is a version of the famous competitive 
exclusion principle (Armstrong and McGehee 1980). 

As pointed out by Nicolis and Prigogine (1977), in the fluid at the onset of 
convection, fluctuations in the form of convection cells appear. These cells or 
rolls can have different sizes, parametrized by the geometric factor a. Solving the 
Lorenz dynamics does not allow us to calculate the size of the convection rolls, 
because a is treated as a constant parameter. But as the monoculture resource- 
consumer ecosystem can be generalized to a polyculture resource-consumers 
ecosystem, it is tempting to perform a translation in order to construct a gener- 
alization of the simplified Lorenz system. This adds a new element in the fluid 
systems: Rolls with different sizes (different a^) will go into competition with 
each other. 

With this generalization, one can can now ask which kind of convection 
cells are the most fittest, which species of rolls will eventually survive. As the 
competitive exclusion principle states, the rolls with the highest fitness fi will 
survive, so the only thing we need to do is to translate the competitive fitness 
measure fi to the fluid system and write it as a function of the parameter . If 
we do the translation with the above dictionary (|10j|16p , we get the fluid fitness 
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for rolls with parameter 



[af + ly-n-^vx 

Note that the geometric factor (|^n|l appears in the fitness. There is a trade-off 
between small and large sizes, and the fitness (|34|) is maximal for rolls with 
parameter a.i — 1/V2, and hence with width h/ai = y/2h. As was first shown 
by Rayleigh (1916) using a totally different line of reasoning, this is also the 
experimentally verified size of the convection rolls at the onset of convection. 
Furthermore, using this value for a together with and the definition of the 
Rayleigh number (llSp . we can calculate the critical Rayleigh number Ruc = 
27ti^/A. This is indeed the correct value for the fiuid system with free-free 
boundary condition^ 



7 Predating fluid motion 

A next natural step to take is describing our ecosystem with the addition of 
predators eating the consumers. This leads us to a more speculative idea: 
Is there a possibility for 'predation' in fluid systems? Let us first study the 
resource-consumer-predator ecosystem 

The dynamical equations for the resource is the same as For the con- 
sumer and the predator the dynamics changes to 

jGc = qcR9cR{CR-C%lK,,)Cc-9PcCcCp-dcCc, (35) 

^Cp = qpcgpcCcCp - dpCp. (36) 
at 

There is now a second critical bifurcation point 

^ + dc dp ^ ^3^^ 

ttfl qcpgcR ap, + KAc 

such that for values A" < A°2 we get the previous solutions (|27ll30p . For values 
higher than this second critical concentration level, the former states become 
unstable and the new stable state has a non-zero predator concentration: 

^* qpcgpcap^" , , . 

Op — 7 ; h — — , (oa) 

qpcgpciap + KAC) + gcpdp J<eq 

C*c = (39) 

qpcgpc 

^* ^ qcRgcBMPcap^^ dc_ ^^Q^j 

^ qpcgpcio!R + KAc) + gcpdp gpc 

Moving to the convective fiuid system, we have to study convection patterns 
that appear beyond a second bifurcation point. As shown above, there is a 
first bifurcation from conduction to straight convection rolls. In the straight 
rolls situation, there was only velocity in the x- and z-directions, leading to a 



^This result is non-trivial, as the final words of appendix |B] point out. 
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non-zero kinetic energy for these two directions. This Ekm.xz was shown to be 
related with the consumer concentration. But for certain systems (depending 
on e.g. the Prandtl number), due to the appearance of a vefocity gradient 
in these roUs, there might be changes in the surface tension leading to a new 
instability at a second critical gradient level. This was experimentally as well 
as numerically shown (Clever and Busse, 1987, Getling 1998). At this second 
bifurcation a new pattern arises, from straight rolls to zig-zag rolls or rolls with 
travelling waves in the direction of its rotation axis (the y-direction). In these 
new patterns, there is also a non-zero velocity component Vy in the y-direction, 
leading to a non-zero kinetic energy Ekin,y 

This allows us to propose a conjecture. The Lorenz system was derived by 
simplifying the Navier-Stokes equations in the Boussinesq approximation. By 
taking the lowest modes in an expansion, and performing an approximation, the 
Lorenz system was derived in order to study straight convection rolls. The wavy 
pattern could not be studied with the Lorenz dynamics. The conjecture states 
that by including another mode, a new variable that describes the motion in 
the y-direction, a new set of dynamical equations can be given (after performing 
some approximations to guarantee that the equations are autonomous), and this 
set of equations can be translated into the dynamics of a resource-consumer- 
predator ecosystem. 

More specifically, the hypothesis that one can make is that the predator 
concentration is proportional with the kinetic energy of Vy. The interpretation 
is that the waves are behaving as predators feeding on the velocity gradient (or 
kinetic energy) of the 'consumer prey' rolls, in a similar way as the consumer 
prey rolls are feeding on the temperature gradient (heat energy). 

We did not prove this conjecture at the level of the dynamical equations, 
but one will only give some (intuitive) arguments. 

First, by looking at the advection term in the heat equation (|B-2p . one can 
see that there is a coupling between temperature and velocity, and it is this 
coupling that was proven to be equivalent with the coupling of consumers with 
the resource in the ecosystem dynamics. Now, by looking at the advection term 
in the Navier-Stokes equation (jB-3|) . one can see that there is indeed a coupling 
between different velocity components, so one might expect that this results in 
an equivalent coupling between predators and consumers. 

Second, our conjecture implies that the predator parameters are related to 
the fluid parameters, in a similar way as in (|10m6p . One might intuitively guess 
that e.g. dp ^ v. As can be seen in ([37|) . a term ^^^^^^ appears. As this is the 
ratio of viscosity over conductivity, this term is proportional with the Prandtl 
number (jl9p . The prefactor is dependent on the geometric factor which now 
includes the wavelength. As shown in e.g. Busse (1978), the second critical 
gradient level increases when the Prandtl number increases. This is consistent 
with the increasing behavior observed in p7p . 

A third test for this 'predator - kinetic energy' hypothesis is performed by 
looking at the thermodynamical level. If our conjecture is correct, the total 
steady state heat transport should be related with the ecosystem metabolism 
in the steady state, as in (PS]) . Using (|38p . the latter can be easily calculated 
and is presented in Fig. El We see that for input concentrations above the 
second bifurcation point, when predation is possible, the stable predator state 
has always a lower ecosystem metabolism rate than the unstable consumer- 
only state. Looking for example at the behavior of the Nusselt number (the 
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dimensionless number which is proportional with the total heat transport) in 
the fluid system (Fig. 6 in Clever and Busse, 1987), we can see that for all 
studied parameters, the heat transport in the wavy roll state is indeed always 
lower than in the straight roll state. 




Figure 6: The total steady state ecosystem metabolism F^jyj in function of 
the external gradient A", for the predator ecosystem. One can clearly see the 
existence of two critical bifurcation points. Beyond the second bifurcation, the 
rate in the predator state (blue) is lower than the consumer state (dashed green) . 

Hopefully, one can rigorously proof this correspondence between ecological 
and fluidal predation. This would allow us to require more analytical expressions 
instead of using numerical simulations (Clever and Busse, 1987). Furthermore, 
if this would be possible, we get a new parameter, the wavelength of the zig- 
zag or wavy pattern which might be related with the parameters qpc: 9 PC 
and dp. Perhaps it is possible to derive the experimentally correct wavelength 
(see Pomeau and Manneville, 1980) again from competition and fitness at the 
predator level, because the competitive exclusion principle also works at this 
level (Smith and Waltman, 1995). 



8 Conclusions and further discussions 

We have seen that one can simplify the dynamical equations of a convective fluid 
system into a set of two ordinary differential equations which look exactly the 
same as a simplified resource-consumer ecosystem. Furthermore, there is not 
only a mathematical correspondence in the structure of the equations, but more 
remarkable, there is also a correspondence between physical interpretations. 
This correspondence was then broadened to include competition and predation. 
With these extensions, we have proven or conjectured more connections. 

• We were able to calculate the correct value of the size of convection cells 
with the help of biological competition and fitness. 

• we have translated quantities, processes, energetic pathways,... from fluid 
systems to ecological systems and vice versa. 
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• and we have conjectured a possible explanation of the decrease in energy 
dissipation in the fluid system after a second bifurcation point as being 
related with the appearance of predative behavior. 
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8/7r4(l + a^) 


decay 


dc 




exchange 


OiR., KAC 




flux 


^RiC^R - 




fitness 


qcRgcR/dc 





Table 2: The correspondence 

Table[2]presents the dictionary of the correspondence which is a quantitative 
extension of the bare essentials given in table [TJ There is always unstructured 
gradient degradation, but above a critical level of the gradient, ordered patterns 
or structures appear: living cells and convection cells. This striking analogy that 
we have found between two systems that are at first sight totally different, can be 
casted in the (more general but often vague) language of dissipative structures 
used by Prigogine and co-workers (Glansdorff and Prigogine, 1971, Nicolis and 
Prigogine, 1977). The structured patterns of the RB system are often putted 
forward as prime examples of dissipative structures. It is believed that life also 
behaves as a dissipative structure. Prigogine and co-workers performed quan- 
titative studies of biological systems, but these were mostly restricted to the 
subcellular level. Schneider and Kay (1994) took the correspondence further 
to the ecosystem level, but their discussion was only qualitative, using often 
vague words. Furthermore, in recent decades a new 'maximum entropy pro- 
duction' (MaxEP) school emerged (Kleidon, 2004; Kleidon and Lorenz 2005; 
Martiouchev and Sclcznev, 2006) where it is believed that complex processes, 
including life, tend to maximize the entropy production. With our work, we ex- 
tended the program initiated by Prigogine, Schneider, Kay and others by study- 
ing quantitatively the (thermodynamic) properties of dissipative structures at 
the ecosystem level. This lead to a more exact formulation of the correspon- 
dence, but also a new feature appeared, something which was not studied by 
Schneider et al.: the appearance of 'predative dissipative structures' after a sec- 
ond bifurcation. Not only is Fig. [6] an extension of Fig. |4] (which was shown 
to be equivalent with Fig. 2a in Schneider and Kay, 1994), it also shows that 
the total gradient degradation (by heat transport or ecosystem metabolism) of 
the consumer-predator state is lower than the corresponding unstable consumer 
state. As the gradient dissipation is proportional with the total entropy pro- 
duction, this might be a criticism on the basic hypotheses of the MaxEP-school. 
For high thermodynamic gradients the energy dissipation of the state with 'sec- 
ond level' predative dissipative structures is lower than the state with only 
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'first level' dissipative structures. The predative dissipative structures make the 
system less efficient in degrading the thermodynamic gradient. 

There are many new questions about the fluid- ecological system analogy. 
We were able to determine the size of convection rolls by translating the Dar- 
winian view of evolution and natural selection to the fluid system. Can this 
be generalized, i.e. is the difficult pattern selection problem in fluid systems 
(Getling, 1998) analogous to the difficult problem of evolution and natural se- 
lection in ecology? More specifically: What about 'predatory' pattern selection 
(e.g. the selection of the wavelength of the wavy rolls)? What about turbulent 
fluid states, longer trophic chains (top-predation), 'fluidal niches and food webs', 
evolution at different time scales, genetic information, velocity correlations,...? 
Up till now, we were not yet able to derive new non-trivial results, because solved 
problems were related with solved ones, and unsolved with unsolved ones. We 
hope that besides the esthetically pleasing results we have found, one is able to 
use the analogy to find new solutions to important problems, both in ecology 
and fluid physics. 

A The RCW ecosystem 

The resource-consumer-waste ecosystem consists of two environmental reser- 
voirs, one for the resource and one for the waste. As an example, we can think 
of a chemotrophic ecosystem with glucose or methane as resource and CO2 as 
waste product. The resource is supplied from the environmental reservoir at a 
fixed concentration using a linear exchange mechanism with rate constant 
an and flux 

FR = aR{Cl-CB). (A-1) 

Cii{t) is the variable resource concentration in the ecosystem. The ecosystem 
metabolism is the total conversion (degradation) of resource into waste. In our 
chemotrophic ecosystem, this conversion is an oxidation proces. Table |3] shows 
the three metabolic transformations that occur within the ecosystem, together 
with the kinetic expressions used. 



Abiotic conversion 




Fac = ^^AC [Cp, ~ 


Biomass synthesis 
and biotic conversion 


R ^ qcrC 

+ {l-qcR)W 


Fbc - gcpCc [Cp, - 


Biomass decay 




Fbd = dc'Cc 



Table 3: Ecosystem transformations 

The abiotic conversion is a chemical reaction with equilibrium constant K^q 
and a constant abiotic conversion rate parameter kac- The latter abiotic con- 
version rate is increased due to a parallel biotic conversion, described by a simple 
linear functional response with parameter gcR- This biotic conversion has two 
parts: a fraction of the resource is used for consumer growth, the other part 
of the resource turns immediately into waste. From a thermodynamic perspec- 
tive, the latter resource turnover is necessary to drive the growth process. This 
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fractioning is described by the yield parameter qcr < 1: this is the growth 
efficiency which denotes the amount of resource required to build up one unit 
of biomass. The third metabolic transformation is the biotic decay (biomass 
turnover), represented by the rate constant dc- 

When the resource is turned into waste, the latter is emitted into the waste 
reservoir from the environment. The latter has a constant waste concentartion 
C^, and the exchange flux can be described as 

Fw =awiC^-Cw)- (A-2) 

Putting the two exchange fluxes and the three metabolic fluxes together, the 
complete dynamics for the resource concentration Cji(t), the consumer biomass 
concentration Cc{t) and the waste concentration Cw{t) now look like 

^^Cr = aR{C°i-CR)^{KAc+9cRCc){CR-Cw/K,g), (A-3) 

j^Cc = qcR9cRiCR-Cw/K,g)Cc-dcCc, (A-4) 

—Cw — awiC^ — Cw) 
at 

+ {kac + (1 - qcR)9cRCc){CR - Cw/Keg) + dcCc- (A-5) 

This is the RCW model. Next, we have to simplify this model to the RC model, 
by assuming aw to be very large. This means that the relaxation time of the 
waste exchange is negligibly small, and we get the condition that Cw ~ Cyi/j 
resulting into ([BllT]). 



B The XYZ Lorenz system 

In this appendix, we will give all approximations and a schematic derivation in 
order to arrive at the Lorenz system for the Rayleigh-Benard convective fluid 
(see Berge and Pomeau, 1984 or Lorenz, 1963). 

In order to present the field equations we will first list the Boussinesq ap- 
proximations (see e.g. Getling, 1998): 

• There are no pressure terms in the energy balance equation. 

• The heat conduction coefficient x ^^J^d the kinetic viscosity v are constants. 

• The local density field p depends on the temperature as p = po{l — a{T — 
To)) with po and Tq the constant reference density and temperature, T the 
local temperature field, and a the constant thermal expansion coefficient. 

• The above dependence of the density on the temperature is taken into 
account only in the gravitational force term in the momentum balance 
equation. At other places in the equations, we will write the density as 
Pa- 

• The fluid is incompressible (except in the thermal expansion term): ^ = 
0, which results in an equality between heat capacities at constant pressure 
and volume: Cp = c^,, or it can be written in terms of the velocity field v 
as: 

V-w = 0,. (B-1) 
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• The local internal energy differential is dU = CpdT. 

With these approximations, the heat transport equation can be derived from 
an energy balance equation, and looks like 

dT 

— = -{v-V)T + x^T. (B-2) 

The first term on the right hand side is the advective heat transport term, and 
the second is the heat conduction term. 

The equation for the velocity field is derived from the momentum balance, 
and results into the Navier-Stokes equation. In the Boussinesq approximation, 
this leads to 

^ = _(^;.V)t;-^-5^1, + i.Az/, (B-3) 
at po po 

with p the pressure field, g the gravitational acceleration and 1^ the unit vector 
in the vertical z-direction. On the right hand side we see respectively the ad- 
vection term, the pressure gradient term, the external gravitational force term 
and the viscous diffusion term. 

As a final step, in order to fully describe our system, we need boundary 
conditions. The boundary condition for the temperature is simply 

T(z = 0) = T°, (B-4) 

Tiz^h) = T". (B-5) 

For the velocity, we have 

v,{z^O) = 0, (B-6) 

v^iz^h) = 0, (B-7) 

because there is no fluid flowing out of the layer. This is not enough, and we need 
another condition on the velocity. We will take free-free boundary conditions to 
make the description of the solutions easier. This gives 

^U=o.H=0. (B-8) 
oz 

In summary, we have five partial differential equations: Three from the three 
velocity components, one from the incompressibility condition and one from the 
temperature. Our five local variables are the velocity, pressure and temperature 
fields. Lorenz made some further assumptions in order to turn these five p.d.e.'s 
into three o.d.e.'s with only three global variables. 

Due to (|B-ip . one can write the velocity field as v — W x tp, with tp the 
streamfunction. We know from experiment that at the onset of convection 
(near the critical gradient), a convection roll pattern will arise (Getling, 1998). 
Suppose that the axis of the rolls are along the horizontal y-direction. Hence, 
there will be no Vy component. The simplest way to obtain this is by assuming 

i>x = ipz ^ 0. 

Next, we want to circumvent the pressure field. This can be done by taking 
the curl of the velocity equation, resulting into: 

_ ^^_^^ + .V^(V^^) + ,«^ (B-9) 
dt ~ dz dx dx dz iV ^/;,j+ga^^. {a M) 
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As a final step, we will expand the temperature and ipy fields in Fourier 
modes, taking the boundary conditions into account, and we will retain only 
three of these modes: 

= X(t)sin(^)sin(^), (B-10) 

T = T^-/3z + y(i)cos(^)sin(^)-Z(t)sin(^), (B-U) 

n n n 

with the width of the convection cell equal to h/a. In Fig. [T]a physical interpre- 
tation is given to the variables X, Y and Z. Plugging these expressions into the 
above partial differential equations (|B-2[) and (jB-9[) . and collecting the factors 
with the same spatial dependence, gives: 

dX _ v{a^ + l)7r^ gaha 

at TT(a'^ + 1) 

-dt - — ^ + — (^-1^) 

dZ _ ^4^2 

As can be seen, the system does not close because there is a cos(27rz//i) term. 
A final approximation consists of taking this cosine equal to one. 

We finally arrive at the Lorenz equations, which we will call the XYZ model. 
Next, we have to simplify this XYZ model to the XZ model, by assuming the 
pseudo steady state condition for Y (i.e. taking dY/dt = 0), resulting into 

We conclude this appendix with an important remark. There are two im- 
portant approximations for the XZ model. The first is the cancelation of the 
cosine factor in (|B-13p . Therefore, solutions of the Lorenz system are not exact 
solutions of the complete fluid system in the Boussinesq approximation. In this 
sense, the result of section [5] is not trivial, because we arrived at the correct 
answer whereas the underlying dynamics does not give exact solutions. 

Our second approximation is the pseudo steady state restriction. This means 
that the steady states of the XZ system are also steady states of the XYZ model 
(but as mentioned above, not necessarily of the complete fluid system). In this 
article we mostly restricted the discussion to the steady states of the XZ model, 
but one should be cautious to use this model to try to find correct transient 
solutions for the XYZ or the complete fiuid systems. As an example, the XYZ 
system has chaotic solutions with unstable steady states, whereas these chaotic 
solutions are absent in the XZ model. The latter has always stable steady states. 
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